Microcanonical fluctuations of the condensate in weakly interacting Bose gases 
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We study fluctuations of the number of Bose condensed atoms in weakly interacting homogeneous 
and trapped gases. For a homogeneous system we apply the particle-number-conserving formulation 
of the Bogoliubov theory and calculate the condensate fluctuations within the canonical and the 
microcanonical ensembles. We demonstrate that, at least in the low-temperature regime, predictions 
of the particle-number-conserving and traditional, nonconserving theory are identical, and lead to 
the anomalous scaling of fluctuations. Furthermore, the microcanonical fluctuations differ from the 
canonical ones by a quantity which scales normally in the number of particles, thus predictions of 
both ensembles are equivalent in the thermodynamic limit. We observe a similar behavior for a 
weakly interacting gas in a harmonic trap. This is in contrast to the trapped, ideal gas, where 
microcanonical and canonical fluctuations are different in the thermodynamic limit. 

PACS numbers: 03.75.Hh, 05.30.Jp 



I. INTRODUCTION 

Experimental realization of atomic Bose-Einstein con- 
densation and degenerated Fermi gases, has stimulated a 
large interest in the physics of ultracold gases [1| . Among 
others, the issue of fluctuations in the number of con- 
densed atoms has been a subject of intensive theoretical 
studies 0]-[l3- ^ is well-known that the standard text- 
book approach based on the grand-canonical ensemble 
predicts unphysically large fluctuations. This pathologi- 
cal behavior is avoided, however, when one employs the 
canonical or the microcanonical ensemble, where the to- 
tal number of particles is fixed. From the experimental 
point of view, trapped atoms are under conditions of al- 
most complete isolation. This favors the use of micro- 
canonical ensemble, at least in the situations when pre- 
dictions of different statistical ensembles are not equiv- 
alent. While in ideal gases the microcanonical fluctua- 
tions have been thoroughly investigated O-jg, the cor- 
responding problem for interacting gases has been only 
studied in the first-order perturbation theory [ll| , or for 
relatively small systems by means of numerical methods 
|12|. On the other hand, the Bogoliubov theory |19| . 
which has proved to be immensely successful to describe 
weakly interacting gases, has been applied to study the 
fluctuations in the canonical ensemble [l(J, Ej, uM , but 
to our best knowledge, it was not used to investigate the 
microcanonical fluctuations. 

The purpose of this paper is to address the problem of 
microcanonical fluctuations in the framework of Bogoli- 
ubov theory. For a homogeneous system we perform our 
analysis within the particle-number-conserving formula- 
tion of the Bogoliubov theory, developed by Girardeau 
and Arnowitt [2l| [35( . In this way we avoid the possible 
errors resulting from the nonconservation of the number 
of particles in the Bogoliubov method. We obtain that 
the microcanonical fluctuations, similarly to the canoni- 
cal ones |lfjt Il4l Il5| , exhibit anomalous scaling with the 
total number of particles. Moreover, fluctuations in the 
microcanonical and the canonical ensemble differ by an 



extensive quantity, and predictions of both ensembles be- 
come equivalent in the thermodynamic limit. For a ho- 
mogeneous system this feature is already observed for an 
ideal gas. Nevertheless, in a harmonic trap, the micro- 
canonical fluctuations of noninteracting condensate re- 
main different from the canonical ones even in the ther- 
modynamic limit. We show here, that inclusion of weak 
interactions described in terms of Bogoliubov theory, re- 
stores the equivalence of microcanonical and canonical 
description, and gives rise to the anomalous scaling of 
fluctuations in both considered ensembles. Noteworthy, 
for a homogeneous system we observe that application of 
the particle-number-conserving method leads to the sim- 
ilar results as obtained from the standard, nonconserving 
approach. Therefore, for a trapped gas we carry out our 
derivations only within the nonconserving theory, which 
we find more convenient in this case. We discuss this is- 
sue in the last section explaining why at low temperatures 
the fluctuations are insensitive to the actual number of 
atoms. 

The anomalous scaling of the fluctuations predicted 
by the Bogoliubov theory is closely related with the 
existence of phonon excitations at low energies. On 
the contrary, application of the first-order perturbation 
theory results in the normal scaling of the fluctuations 
(llL Il3 . fl5J | . While the perturbative treatment is more ap- 
propriate for finite size systems with sufficiently small in- 
teractions, in the thermodynamic limit the system should 
be rather described within the Bogoliubov theory, since 
the presence of arbitrary small interactions leads to the 
phonon excitations at low energies in the thermodynamic 
limit. We note that there is some controversy on using 
the Bogoliubov method to calculate the quantities higher 
than the second order in the field operators [I3. An ex- 
cellent agreement of the Bogoliubov theory with the ex- 
act statistics calculated for the one dimensional trapped 
gas |lq | provides, however, strong arguments in favor of 
the Bogoliubov approach. The ultimate test of the the- 
oretical predictions should be provided by experiments, 
either through the measurement of the second-order cor- 



relation functions, or by detecting the statistics of scat- 
tered photons 20] . 

The paper is organized as follows. In section [HI for 
the sake of completeness, we analyze the fluctuations 
of ideal Bose gases confined in a three dimensional box 
and a harmonic trap. The behavior of fluctuations in 
an interacting gas confined in a three dimensional box is 
studied in section ITTT1 In particular, in section UlI Al we 
briefly describe the Girardeau- Arnowitt particle-number- 
conserving formalism. Section 1111 Bl is devoted to the 
analysis of fluctuations in the canonical ensemble, while 
the microcanonical fluctuations are calculated in section 
IIII CI Section Hvl presents the results for a weakly inter- 
acting trapped gas. The canonical fluctuations are con- 
sidered in section IIV Al whereas section IIVBI describes 
the results for the microcanonical ensemble. We end in 
section presenting some conclusions. Finally, four ap- 
pendices give some technical details. 



II. IDEAL GAS 

In this section we analyze fluctuations of the num- 
ber of condensed atoms in an ideal gas. First, we con- 
sider homogeneous system with N noninteracting bosons 
confined in a box of the size L with periodic boundary 
conditions. Below the critical temperature of the Bose- 
Einstein condensation, the fluctuations of the number of 
condensed atoms can be calculated with the help of the 
Maxwell's Demon ensemble (MDE) 4]. In this approach 
the condensate serves as an infinite reservoir of particles 
for the subsystem of excited states. Such treatment is 
justified as long as the probability of states with all the 
atoms excited is negligible, hence it is not valid in the 
vicinity of the critical temperature. The mean number 
of excited atoms (N e ) and its fluctuations (8 2 N e ) in the 
canonical ensemble, can be calculated from the grand 
canonical partition function S e (z, (3) of the excited sub- 
system Q 



dz 



(N e ) CN = z— lnS e (z,/3) 
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In this formula e v denote the single-particle energy of 
the level u, f3 = 1/fcgT is the inverse temperature and 
z is the fugacity. In the canonical ensemble, the to- 
tal number of particles A" is conserved, and obviously 
<^o) CN = N - <iV e ) CN , and (S 2 N ) CN = (S 2 N e ) CN , 
where (Ao) CN and (<5 2 Af ) CN denote the mean number 
of condensed atoms and its fluctuations, respectively. In 
the thermodynamic limit, summation on the r.h.s. of 
Eqs. QJ and J2J, can be replaced by integration. How- 
ever, this is not always possible, in particular, in the 



case of three dimensional box, where resulting integral 
for the condensate fluctuations is infrared divergent. To 
avoid these difficulties, in our study the thermodynamic 
quantities are evaluated from the integral representa- 
tion involving spectral Zeta function 0, |8j. Employ- 
ing integral representation of the exponential function 
e~ f = (2iri)^ 1 J _. dz t~ z T(z), one can express (N e ) 
and (<5 2 A r e ) as contour integrals of Mellin-Barnes type 
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— Z(f3,z)C(z)T(z), (3) 
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(S 2 N e ) c ^ = / — Z(p,z)az-l)T(z), (4) 
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where £(z) denotes Riemann Zeta function, and Z(/3, z) is 
the spectral Zeta function: Z((3,z) — 2~2v^o(@£v)~ z ■ The 
contour of integration lies to the right of all the poles of 
the integrand. Now, the problem of calculating fluctua- 
tions has been reduced to determination of the poles of 
Z(/3,z). For the considered case of a three dimensional 
box with periodic boundary conditions, the spectral Zeta 
function reads 



Z(J3,z) 



Ce(z) 
G9A)« 



(5) 



where Ce(z) denotes the Zeta function of Epstein [24 



Ce(z) = £ 



{n% + n* + nlY ! 



Rez > -, (6) 



and A = 2tt 2 H 2 /mL 2 is the energy of the first excited 
state. Function (e(z) has merely single pole at z = 3/2 
equal to 2tt. The straightforward calculation of the con- 
tour integrals J2J and JIJ with Zeta function j£|, leads 
to 



(S 2 N e ) 



CN 



CN 



TT 3 / 2 C(3/2)t^ 2 + C E (l)t, (7) 

C £ (2)i 2 +^ 3 / 2 C(l/2)t 3 / 2 , (8) 



where t — ksT/A. In derivation of Eqs. Q and (JSJ, we 
have included leading and next to the leading rightmost 
pole of the integrand. Thus, the result accounts also for 
correction due to the finite size of the system. As was 
first pointed out in |25|. the canonical fluctuations in a 
box are anomalous. This can be observed by rewriting 
Eq. JSJl in the form containing explicit dependence on 
the number of particles 
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CN 
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where f = T/T c , and k B T c /A = n-^N/Ci^)) 2 ^ is the 
critical temperature. 

Now we turn to the microcanonical ensemble. In or- 
der to calculate fluctuations of the number of condensed 



atoms, we employ the following relation 



W M C = W 



[(SNeSE)^] 



CN 



(S 2 E) 
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CN 



which links microcanonical fluctuations with the quanti- 
ties calculated in the canonical ensemble: particle-energy 
correlation {6N e 5E) C N = ((N e -{N e ))(E- (£}))cn, and 
fluctuations of the system's energy (5 2 E)qn = ((E — 
(i?)) 2 )cN- The former can be determined from the mean 
number of excited particles (N e ) 



(SNJE) 



fdN e 



CN 



V dp 
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In a similar manner, one can calculate the fluctuations of 
energy from its expectation value 



(6*E) 



CN 



dE 
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while the mean energy (E) can be expressed as a contour 
integral containing spectral Zeta function Z(/3, z) 
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After straightforward calculations, we arrive at the fol- 
lowing result for the microcanonical fluctuations in a box 
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By comparison with the canonical fluctuations (jSJ), we 
see that the leading-order term, which gives rise to the 
anomalous scaling, is the same. The difference between 
fluctuations in the considered ensembles appears in the 
finite size correction term, thus the microcanonical and 
canonical fluctuations become equal in the thermody- 
namic limit. 

Let us turn now to the case of a trapped gas. For 
simplicity we assume that the atoms are confined in a 
spherically symmetric harmonic potential, however gen- 
eralization of the results for anisotropic traps is straight- 
forward. The spectral Zeta function for the harmonic 
trap of frequency u is given by |8j 

Z(/3, z) = r z (*C(* - 2) + f CO* - 1) + C(*)) - (15) 

with r = ksT '/(fko). Substituting this function into 
Eq. Q and taking into account leading and next to the 
leading rightmost poles, we obtain the following result 
for the fluctuations in the canonical ensemble 

(S 2 N e ) CN = C(2)r 3 + r 2 (f lnr + C(2) + | + | 7 ) , (16) 

where 7 ~ 0.577 denotes Euler's constant. The logarith- 
mic correction in l|l()|) originates from a double pole of 
the integrand at z = 2. In the thermodynamic limit: 



FIG. 1: Fluctuations of the number of excited particles in 
the microcanonical and canonical ensemble for the system of 
N = 1000 noninteracting bosons confined in a box with pe- 
riodic boundary conditions (upper panel) and in a harmonic 
trap (lower panel). Analytical results (solid line) derived in 
MDE approximation are compared to the exact numerical 
data (dashed line). The temperature is scaled in the units 
of the critical temperature Tc- 



N — > 00, u> — ► and Nu> 3 —const, the leading order term 
in Eq. I|ltjfl is proportional to N and the fluctuations ex- 
hibit normal behavior. In the microcanonical ensemble 
we employ the identity (|10fl , which leads to the following 
result for the fluctuations 
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We observe that the leading order term, which gives rise 
to the normal scaling of fluctuation in the thermody- 
namic limit, is different than in the canonical ensemble. 
Hence the microcanonical and the canonical fluctuations 
in the harmonic trap remain different even in the ther- 
modynamic limit, in contrast to their behavior in a ho- 
mogeneous system. 



Fig.^presents fluctuations of the number of condensed 
atoms in the canonical and the microcanonical ensembles, 
calculated for the system of N = 1000 atoms confined in 
a three-dimensional box (upper panel) and a three dimen- 
sional harmonic trap (lower panel). It compares the an- 
alytical predictions of Eqs. JBJ), |Q (box) and Eqs. (fit))) . 
()17J) (harmonic trap) with the exact numerical results cal- 
culated with the help of the recurrence algorithms for the 
canonical and microcanonical partition functions. We see 
that analytic results agree very well with the numerical 
data, apart from the region close to the critical tempera- 
ture, and above it, where the MDE approximation ceases 
to be valid. 



III. WEAKLY INTERACTING 
HOMOGENEOUS GAS 



Girardeau- Arnowitt particle- number-conserving 
formalism 



In this section we present a brief description of the 
particle-number-conserving version of the Bogoliubov 
theory, formulated by Girardeau and Arnowitt (GA) [2l| . 
In the following, we restrict our analysis to the regime of 
temperatures much lower than the critical temperature. 
This allows us to make the following assumptions: (i) the 
mean number of condensed atoms is much larger that its 
fluctuations (Nq) ^> (6 No), thus we can we can exclude 
the possibility of exciting the state where the conden- 
sate is totally depleted, (ii) we neglect the influence of a 
finite temperature on the excitation spectrum, through 
the thermal depletion of the condensate, which is pre- 
dicted in the Popov theory 26] . In addition we assume 
dilute-gas regime na 3 <C 1, where n is the density of 
atoms and a is the s-wave scattering length, which al- 
lows to neglect the influence of the quantum depletion of 
the condensate on the Bogoliubov excitation spectrum. 
Let us introduce unitary operators /3 = (N Q + l) _1 / 2 a , 

$1 = aJ(#o + l)~ 1/a , /3qA> = AjAS = 1. where a k, 4 are > 
respectively, the annihilation and creation operators of 



the particles with momentum hk, and N Q — a do- The 
particle-number-conserving ground state of GA theory is 



with the unitary operator U 



(19) 



U = exp 
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Coefficients -0k are real, even function of k. In the GA 
theory they are determined by minimizing the expecta- 
tion value of the full Hamiltonian in variational ground 
state <|19|) . The inclusion of the total Hamiltonian in 



the variational calculus leads, however, to the excita- 
tion spectrum which possess a gap 2lj, which is not 
appropriate for studying the system in the thermody- 
namic limit. Gaplcss Bogoliubov spectrum is recovered, 
when one minimizes the expectation value of the particle- 
number-conserving Bogoliubov Hamiltonian 
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where g = AirT^a/m, V = L 3 , and N is the opera- 
tor of the total number of particles. This approach as- 
sumes omitting the Hartree-Fock and pair-pair interac- 
tion terms in the part of the total Hamiltonian, which 
gives nonzero contribution in the variational ground 
state. These terms, however, are negligible in the consid- 
ered range of temperatures. The minimization procedure 
yields J2J 



tanh tpk 



£k 



gn - £k 



gn 



(22) 



where e k = h 2 k 2 /2m, n = N/V is the atomic density, 
and e^ = \J e\ + 2gne\^ are the energy levels of the Bo- 
goliubov spectrum. 

The excited states in GA theory are constructed in 
basically similar manner to the variational ground state 
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and N = N — ^w n k- Here, {n k } is the set of integer 
numbers, where n k represents the number of elementary 
excitations with energy e k . 

B. Fluctuations in the canonical ensemble 

In an interacting gas, canonical fluctuations of the 
number of condensed atoms can be calculated with the 
help of the following generating function. 



X(z,0)=Tr 



{/«e-^} 



(25) 



where N e — X^ k =^o ^iL^k is the operator of the number 
of noncondensed particles, and the trace has to be taken 
over all eigenstates of the system with the total number 
of particles equal to N. For z = \, generating function 
X(z,/3) is equal, by definition, to the canonical parti- 
tion function. The fluctuations of the number of noncon- 
densed atoms (5 2 N e ) are given by 



<^> CN =*J^lnX(*,/3) 



(26) 



Obviously, the total number of particles is conserved, and 
the fluctuations in the number of condensed and noncon- 



densed atoms are equal: (S 2 N ) CN = (<5 2 A e ) CN . We cal- tiplied by the factor (u 2 . 



In comparison to the ideal gas result J5J), thermal fluctu- 
ations of the number of elementary excitations, are mul- 
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culate X(z, (3) by performing trace in the basis of eigen- 
states | ^-[Ti k }) °f GA theory 

= y E(*M\*°~ 1 *'*'- f *~ 1 ***\*i»*})> 

{»k} 

(27) 

where summation runs over sets {n^} of the occupation 
numbers of elementary excitations. In principle, accord- 
ing to the definition of 4 , { nk i.), the number of elemen- 
tary excitation X^k^o nk cann °t exceed the total number 
of particles N. However, well below the critical tempera- 
ture, the probability of exciting such a state is negligible, 
and we can consider the summation in Eq. 1|27[) as un ~ 
constrained. Similar approximation is assumed in MDE, 
where we treat the condensate as a reservoir of particles 
for excited subsystem, and as a consequence we consider 
the partitions of excited atoms in the canonical ensemble 
without the particle number constrain. The transformed 
operators C/ _1 7V e [7 and U^ 1 HbU can be found with the 
help of the following identity 



U 1 a\JJ = Wkflk 
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Second term 

on r.h.s. of Eq. (|3 II) . describes the quantum fluctuations 
at zero temperature S 2 N q = 2 J^wg u k v k- We note that 
a similar expression can be derived in the framework of 
the standard, particle-number- nonconserving theory |l(1 | . 
We comment on this issue in the last section of the pa- 
per. Now, we rewrite (5 2 N e } in terms of contour integral 
involving spectral Zeta function. Calculations similar to 
the derivation of J2J and J3J, yield the following integral 
representations 



(S 2 N e ) CN = / —Z(t,z,a)C(z-l)r(z) 



+ S 2 N q , 
with spectral Zeta function 



(32) 



Z(t,z,a) = t 2^ =j+2-, (33) 



where a = gn/A = 2Na/nL, and n = (n x ,n y ,n z ) is 
a vector composed of integer numbers. Spectral Zeta 
function Zi(t,z,a) can be expressed as a sum of two 
simpler functions 



where Wk = cosh V>k and «k = — sinh ipk ar e the usual am- 
plitudes of the Bogoliubov transformation. After some 
straightforward algebra we obtain 



k/0 k#0 

+ Yj ((/3o) _2 aka-k + H.c. 
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tJ- l H B U ~ £ + ^£k4«k, 
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where Eq is the ground state energy of the system. The 
transformation U ■ U does not diagonalize the Hamil- 
tonian Hb exactly, and in the derivation of l|3UI) we have 
neglected terms describing interactions between excited 
atoms. This approximation is consistent with neglection 
of the Hartree-Fock and pair-pair scattering terms in the 
total Hamiltonian, and is well justified in the considered 
range of temperatures. 

Employing Eqs. ©, ©, JH and (HI we obtain 
the following result for the fluctuations in the canonical 
ensemble 
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Z(t, z, a) = t z (9(z, a) + 2a 2 9(z + 2, a)) , 



with 



(z,a) = Y 
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In the thermodynamic limit: A— >oo,V^— >oo,n = 
const, and a ~ n 1 / 3 ]^ 2 / 3 — > oo. Thus, for sufficiently 
large system we can consider the regime a»l, and in 
the subsequent derivation we will approximate function 
6(z, a), by its asymptotic expansion for large parameter 
a. The analytic properties of the function 9{z, a) are 
studied in appendix ^] while its asymptotic expansion is 
developed in appendix IBl In the former it is shown, that 
the poles of 8(z, a) are located at z = | — m, m = 0, 1, . . ., 
with residues equal to 27r((-2a) m /m!)r(| + f )/T(| - 
-?). Hence the residues of Z(t, z, a) are proportional to 
(a/t) m = (gn/kBT) m , and inclusion in the calculation 
only the leading rightmost poles, would yield the result 
which is valid only in the regime of temperatures ksT 3> 
gn In order to obtain result, which holds also for low 
temperatures: kgT < gn, in the calculation of contour 
integral l|32l) , we have to sum over all the poles which lie 
to the left of the contour of integration. This summation 
can be done analytically, the details of the derivation are 
presented in appendix [CJ The result for the canonical 



fluctuations reads 
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where 
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and £f/(s,a) is Zeta function of Hurwitz: ^h(s,ci) = 
Ylk=o(k + a ) _S ; defined by the series for Re s > 1 . The 
leading order term: i^(2)i 2 is proportional to iV 4 / 3 in 
the thermodynamic limit, and the canonical fluctuations 
in the interacting gas are anomalous, similarly to the 
ideal gas. The prefactor of the anomalous term is ex- 
actly twice smaller than for the ideal gas, which can be 
attributed to the coupling of k and — k modes in the 
Bogoliubov Hamiltonian [l4J]. We note that our result 
agree in the leading order with the result of Ref. jlCj . 
calculated in the framework of the standard, particle- 
number-nonconserving theory. 

It is interesting to investigate the behavior of fluctu- 
ations in the two limits: fc^T 3> gn and fc^T <C gn. 
For the assumed range of temperatures, the former can 
be realized in the system with sufficiently small interac- 
tions. The latter corresponds to the regime where the 
thermal excitations occur mainly in the phonon part of 
the energy spectrum. In the limit fcgT 3> gn, we perform 
Taylor expansion of functions Ch(s, 1 + n) and si(j]) for 
small 7/ = a/2-jrt, and obtain 
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In the regime k B T <C gn, we perform asymptotic ex- 
pansion of the functions £h (s, 1 + 77) and Si(rj) for large 
n = a/2-rrt J22J, and obtain the following result 
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FIG. 2: Fluctuations of the number of noncondensed atoms 
in the canonical ensemble for interacting Bose gas confined 
in a box with periodic boundary conditions. Plots are shown 
for a = 2.4 and a — 10 (inset), where a — 2Na/nL. The 
analytical result of Eq. 136H is compared with its low and high 
temperature expansions: Eqs. I39H and (I38L respectively. 
The temperature is scaled in the units of A = 2ir 2 h 2 /mL 2 . 



Fig. [3 presents the fluctuations of the number of non- 
condensed atoms in the canonical ensemble calculated 
for a — 2.4. The inset shows the same quantity for the 
system with stronger interactions: a = 10. The high- 
est temperature presented in the plot, corresponds to 
T/T c « 0.5 for N = 1000. Fig. H compares the analyti- 
cal result of Eq. I|36|) , with its low and high temperature 
expansions: l|59")l and (|35)l . respectively. We observe that 
the low temperature expansion correctly describes the 
fluctuations in the whole regime of considered tempera- 
tures. This reflects the fact that the main contribution 
to the fluctuations, comes from the phonon part of the 
spectrum. On the other hand, the high temperature ex- 
pansion is valid only for temperatures t > a. 



C. Fluctuations in the microcanonical ensemble 
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CN 



\W2)t 2 



TrV2at+— (2a) J 
4 



(39) 



This result can be also derived, by assuming phonon dis- 
persion relation for the energy spectrum: e^ — ck, with 
the sound velocity c = Wgn/rn . In this case spectral 
Zeta function takes form 



Z p h(i, z, a) 



^^2 



1 



(40) 



The last equation can be obtained directly from Eq. (|33|l , 
by approximating each term of the series by its asymp- 
totic behavior for a — > 00. 



We now turn to the the microcanonical ensemble. 
Fluctuations of the number of noncondensed atoms can 
be calculated from the following generating function 



Y(z,E)=Tr{z*°5(E-H B )} 



(41) 



where delta operator 5(E — Hb) selects the states with 
the total energy equal to E, and the trace is taken over 
all the eigenstates with a fixed total number of particles. 
Taking into account the definition of Y(z,E), the mean 
number and fluctuations in the microcanonical ensemble 



can be evaluated from 



(N e 



MC 



(S 2 N e ) 



MC 



z—\nY(z,E) 
dz 



z—z—lnY(z,E) 
oz oz 



(42) 
(43) 



According to Eqs. 142|) and (|43[1 . the fluctuations of the 
number of noncondcnsed particles can be expressed as 



(5 2 N e ) 



MC 



d(N, 



e/MC 



dz 



(44) 



where derivative of the mean number of noncondensed 
particles should be taken before setting z = 1 in Eq. 
(|42[l for (N e ) MC . In a similar manner, one can represent 
the fluctuations of the number of noncondensed particles 
in the canonical ensemble 



<<^) CN = Z ( 



d(N, 



e/CN 



dz 



(45) 



Now we are ready to establish the relation between fluc- 
tuations in the microcanonical and canonical ensemble. 
In the following we show that identity (fTU|) , derived pre- 
viously for an ideal gas, holds also in the case of an in- 
teracting gas. We basically repeat the steps of derivation 
for the ideal gas. First, we assume that the mean num- 
ber of noncondensed particles in both considered ensem- 
bles coincide. This must be true in the thermodynamic 
limit, where we expect that both ensembles are equiva- 
lent with respect to the basic thermodynamic quantities. 
According to this assumption, we denote the mean num- 
ber of noncondensed particles by N e , not distinguishing 
between microcanonical and canonical ensembles. We 
choose energy E and fugacity z as a pair of state vari- 
ables and take the total differential of N e 



dN P = 



fdN e 



V dz 



dz 



fdN e 



V dE 



dE. 



From the total differential dN e it is straightforward to 
derive the following relation 



[8N e 



\ dz 



fdN e 



\ dz 



fdN e 



\dE 



f I <«) 



According to Eq. l(I5|. l.h.s. of Eq. (g7|) at z = 1 
represents the canonical fluctuations of the number of 



noncondensed atoms. The r.h.s. of 147(1 can be rewritten 
in the following way 



8N e 



dz ) T \ dz / E 



Now, if we observe that 



dN e \ «(a*^®), 



^ 2 (If) 



(48) 



knT 2 



dNe 
dT 



z=l 



dE 
~dz~ 



z=l 



and 



(8N e SE) CN , 
(49) 



k R T 2 



dE 
~d~T 



(6>E) 



CN 



(50) 



we arrive at the desired relation stated by Eq. i|10fl . 
Thus, we have expressed the microcanonical fluctua- 
tions in terms of the quantities calculated in the canon- 
ical ensemble. Moreover, the particle-energy correlation 
(dN e SE) CN and fluctuations of the total energy (<5 2 15) 
that enter this relation, according to Eqs. H49|) and 
(|50|1 can be determined directly from the mean number 
of noncondensed particles N e and the mean energy E. 
We calculate the latter quantities in a similar manner to 
the canonical fluctuations, representing them in terms of 
contour integrals containing spectral Zeta function, and 
summing the contributions from all the poles of the in- 
tegrand. This yields 



(N e ) = -47rMOf (i ! + »?)+ (|Ce(1) - ^ 2 ^) 



t, 

(51) 



and 



(46) (E/A) = - [rf( H (^ l + r,)- 2r& H (-\, 1 + rj) 



+ Ch(-1 1 + r,)] SttM + ^ - f (2a) « - t, 

(52) 

where r) = a/27rf, and Eq is the ground-state energy. 
Taking into account the leading order terms in Eqs. i|51|) 
and H52|) . after some algebra we obtain the following re- 
sult for the difference between the canonical and the mi- 
crocanonical fluctuations 



(5 2 N e ) CN -(S 2 N e ) uc = -4. 3 / 2 



(Kh (-i l + v) - vCh (§, 1 + v) + Vv) 2 t 3/2 



VOf (f , 1 + V) + ^V 2 Ch (i 1+v) - 36t?Ch (-§, 1 + v) + 20Ch (-§, 1 + v) 

I 



(53) 



Since £ 3 / 2 = (fc^T/A) 3 / 2 ~ 7VT 3 / 2 , and r\ remains con- fluctuations scales normally with the number of parti- 
stant in the thermodynamic limit, the difference between 



cles. Thus, the microcanonical fluctuations, similarly to 
the canonical ones, exhibit anomalous behavior. Fur- 
thermore, the relative difference between (<5 2 iVe) CN and 
(5 2 Ne) MC tends to zero in the thermodynamic limit, and 
fluctuations in both ensembles become equal for large N. 
Let us investigate the behavior of microcanonical fluc- 
tuations in the two limits: kgT 3> gn and fcgT C gn. 
To this end, we can expand Eq. Q53|l for large and small 
values of n. In the regime a/2irt <C 1, the leading behav- 
ior of fluctuations in the canonical and microcanonical 
ensemble is governed by 






Ce(2) 

Ce(2) 



2q 

2 



t 2 + C(i)(7rf)i, (54) 



.3/2 C(a) 5^(5) 



C(|) 



2j^3/2 



(55) 
(56) 



We notice that in this limit, the difference between en- 
sembles is given in the leading order by the ideal gas 
result (cf Eqs. (JSJ and (|14H h In the opposite limit: 
a/2irt 3> 1, we obtain 



(S 2 N e ) 
(S 2 N e ) 



CN 



MC 



\C E {2)t 2 -nV2^t+\^ 2 {2a)^ 1 (57) 

iC£(2)t 2 -i|vrV2^t+^ 2 (2a)i. (58) 



In this regime, the microcanonical fluctuations differ 
from the canonical ones by the prefactor of the finite 
size correction term. 

In order to verify predictions of the analytic results, we 
performed numerical calculations in the microcanonical 
ensemble. The microcanonical partition functions were 
computed with the help of recurrence algorithms 6, 28] . 
From the partition functions, we calculate the statistics 
of elementary excitations, and finally the fluctuations of 
the number of noncondensed particles from 



(^>MC = 



E 

k,q#0 



[(t 



2 2 \ 

i q + vj 



x ((n k n q ) M c - (nw) Mc(nq)Mc)] 
- 4 E u k u k ((%n-k>MC + (%)mc + \) 

k#0 



(59) 



where h^ = a^a^. This relation can be derived starting 
from the definition of Y(z, E), and performing differenti- 
ation according to Eq. i|43|l . The mean number (nk)MC 
of elementary excitations with energy e^ in the micro- 
canonical ensemble is given by 



{nk)MG = 



£{«„} < $ W}I WE U-'HbU) |$ W} ) 



Y(1,E) 



(60) 

Analogous expression holds for the product of the occu- 
pations numbers (nk'iq)MC I n order to express the final 




FIG. 3: Fluctuations of the number of noncondensed atoms 
in the canonical and microcanonical ensemble for a three di- 
mensional box with periodic boundary conditions. Plots are 
shown for a — 2.4 and a — 10 (inset), where a = 2Na/iiL. 
Depicted are: the analytical result of Eqs. 1361 and I53H 
for the canonical and microcanonical fluctuations, respec- 
tively, the canonical fluctuations calculated numerically from 
Eq. 131H . and the microcanonical fluctuations evaluated nu- 
merically with the help of recurrence algorithms. The tem- 
perature is scaled in the units of A = 2-k 2 h 2 /mL 2 . 



result in terms of the temperature, we calculate the mi- 
crocanonical temperature Tmc from the entropy of the 
system: 1/Tmc = dS(N 7 E)/dE. In our case, the en- 
tropy is given by S(N 7 E) = fee logY(l, E), where the 
r.h.s. depends on the total number of particles only 
through the energy spectrum |36j. 

Fig. |3 presents fluctuations of the number of noncon- 
densed particles in the canonical and microcanonical en- 
semble, calculated for a = 2.4 and a = 10 (inset). The 
highest temperature presented in the plot corresponds 
to T/T c ~ 0.5 for N = 1000 atoms. The figure shows 
the analytic results of Eqs. H^bj) and l|53|l for the canon- 
ical and microcanonical fluctuations, respectively. They 
are compared with the canonical fluctuations calculated 
numerically from Eq. (|31|l and with the microcanoni- 
cal fluctuations evaluated from the recurrence relations 
for the microcanonical partition functions. The analyt- 
ical predictions for both ensembles agree very well with 
the numerical calculations. We note that, the missing 
low temperature part of the numerical curve for the mi- 
crocanonical fluctuations is due to large uncertainty in 
determination of the microcanonical temperature at low 
energies of the system. On the other hand, we see that 
in the system of a moderate size (N ~ 1000) the canon- 
ical and microcanonical fluctuations can be clearly dis- 
tinguished. 



IV. WEAKLY INTERACTING TRAPPED GAS 
A. Fluctuations in the canonical ensemble 

In this section we calculate the canonical fluctuations 
of a weakly interacting gas confined in a spherically sym- 
metric harmonic trap of frequency w. For a homoge- 
neous system we have demonstrated that the particle- 
number-conserving approach leads to the same value of 
fluctuations as the traditional, nonconserving Bogoliubov 
method. Since the trapped gas is more complicated to 
treat analytically in the following we apply the standard 
Bogoliubov theory generalized for an inhomogeneous sys- 
tem, limiting our analysis to the leading order behavior in 
the thermodynamic limit. Similarly to the homogeneous 
case we restrict our considerations to the temperatures 
much smaller than the critical temperature, neglecting 
the influence of the thermal depletion on the excitation 
spectrum. In addition, we assume the dilute-gas regime, 
which allows to neglect the effects of the quantum deple- 
tion. Furthermore, we assume that the number of atoms 
is sufficiently large: Na/a^ a ^> 1, where o,ho = \/h/(muj) 
is the harmonic oscillator length, and the condensate can 
be described in the Thomas-Fermi (TF) regime. This 
is particularly convenient since in this regime the an- 
alytic solutions of the Bogoliubov-de Gennes equations 
are known |29j . 

In the Bogoliubov theory the field operator is decom- 
posed into the condensate wave function tyg and the non- 
condensed part &': \P(r) = x I'o( r ) + vfr(r)'. The approxi- 
mate Hamiltonian with neglected third- and fourth-order 
terms in 'J/' is diagonalized by the Bogoliubov transfor- 
mation 



$'(r) = £(u„(r)S„-«:(r)at 



(61) 



where &„, b J, are the creation and annihilation opera- 
tors of the elementary excitations, while the amplitudes 
u v (r), v v (y) are solutions of the Bogoliubov-de Gennes 
equations. In the TF regime, the condensate density has 
typical "inverted parabola" shape with the condensate 
wave function ^o(r) = ^///(l — r 2 /R 2 )/g, which van- 
ishes outside the condensate radius R = ^2^/(muj 2 ). 
The chemical potential is determined by the normaliza- 
tion of ^o(r) and in the TF approximation is given by 
|U = (15Na/aY LO ) 2 ^ 5 huj/2. The low-energy collective exci- 
tations of the condensate in the TF regime have energies 
El HI 



E n i = TiLOEni = hui{2n + 2nl + 3n + I) 



1/2 



(62) 



which is valid for E„i <C [i- On the other hand, the 
amplitudes u v and v v of those modes are described by 



Jnlm± — 



V4n + 21 + 3 
rJT 2 



i-y 



21 ±1/2 



£,£nl 



X y<P«+W»(l-2y*)Y lT , 



(63) 



where /,,± = u v ± v v , y = r/R, £ = 7kj/(2//), 
Pn ' (x) denote the Jacobi polynomial and Y[ m (9, (j>) 
is the spherical harmonic. In the above formula n = 
0, 1, 2, . . . is the radial quantum number, I and m are the 
quantum numbers of the angular part of the wave func- 
tion, and I =/= when n = 0. Since we consider only the 
lowest collective modes with energies E n i <C /i, in further 
calculations we can neglect the amplitude /„_, and use 
the following approximation: u v s=s v v ~ f v +/2. 

We calculate the fluctuations substituting the oper- 
ator of the number of noncondensed particles N e — 
f d 3 r&(r)'V (r)' into (5 2 N e ) = (N 2 )-(N e ) 2 . This yields 



(S 2 N e ) 



CN 



uX 



J2[ d 3 rut(r)u x (r)) h v h x (64) 



where we retain only the leading order term, and h v — 
(eP e " — 1) _1 . In derivation of Eq. (|64(l we assume that 
elementary excitations are populated according to the 
grand-canonical statistics: (ni/hx)cN — ("i/)cn(Aa)cn 
for v =+- A, (hl)cN = 2(n„)cN + (^)cn, with h v = b\b v . 
This is equivalent with using the MDE for an ideal gas, 
and results from assumption that at sufficiently low tem- 
peratures the number of elementary excitations can be 
considered as unconstrained. Furthermore, in analogy 
to the homogeneous gas we expect that the main con- 
tribution to the fluctuations comes from the low-energy 
phononlike modes, and in the thermodynamic limit we 
can approximate h n i by {fie n i)~ x . The overlap integral 
J d 3 ru*(r)u>(r) is calculated in Appendix ID1 Substitut- 
ing the result of integration into Eq. (|64|l . we arrive at 
the final result for the canonical fluctuations [lj| 



(S 2 N e ) 



CN 



M (S)' 



k B T 



(65) 



with 



-* = E 



2{n+l) 2 (n + l+%) 2 




el l ^n + l + ^) 2 (2n + l + l) 2 \ 



(66) 



In the above equation n, I = 0,1,2,... and I 7^ for 
n = 0. The numerical calculation of Eq. (|66|) yields 
A ~ 0.56. We observe that in the thermodynamic limit 
(N — > 00, ijj — > 0, and Nui 3 — const) the fluctuations 
exhibit anomalous scaling, similarly to the behavior in 
a homogeneous system. Moreover, the T 2 temperature 
dependence of fluctuations becomes universal, both for 
the homogeneous and trapped gases. 



B. Fluctuations in the microcanonical ensemble 

To calculate the microcanonical fluctuations in a trap- 
ped gas we employ the thermodynamic relation l|10|) . In 
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section IIII CI we have shown that this identity applies 
also to an interacting gas. The proof is quite general, 
and in the case of a trapped gas the only difference with 
respect to the homogeneous system is the form of the 
operators N e and Hb, that enter the definition of X(z, (3) 
and Y(z,E). 

According to Eqs. Ijllfl and i|12|l . the particle-energy 
correlation and the fluctuations of the system's energy, 
can be directly calculated from the mean number of non- 
condensed particles and the mean energy. The former 
quantity can be evaluated as the expectation value of 
the operator N e = Jd 3 r\£ r '>'(r)'\I>(r)', where in place of 

^(r)' we substitute decomposition (|MT)l . For ( « 1 we 
use the approximation u v w v v ~ f v +/2, obtaining the 
following result 



(N e 



= 2^/i„ /d 3 i 

1/ J 



K(r)r+iV qd) 



(67) 



where A q d stands for the quantum depletion. Since N q <\ 
does not depend on the temperature, its exact value is 
not important for the subsequent derivation. Now, we 
utilize the results of Appendix iDl for the overlap integral 
between the amplitudes u u , and replace the summation 
by integration, which is valid in the thermodynamic limit. 
This yields H3 



(N e 



^l 



k B T 



N, 



qd- 



(68) 



A similar method cannot be applied, however, to cal- 
culate the mean energy of a weakly interacting trapped 
gas. In this case main contribution comes from the 
boundary of the condensate, where the excitations have 
single-particle character. Derivation based on the semi- 
classical formulation of the Hartree-Fock-Bogoliubov the- 
ory, which includes the effects of collective modes inside, 
and the single-particle excitations outside the conden- 
sate, yields the following law J31J 



(E) = -nN + 



20cm , mV 2 fk _ 



-hjj i , 

/it \ TUjJ I 



,'T 



two 



7/2 



(69) 



Now, with the help of results l|68|) and l|69|l we calculate 
the difference between the fluctuations in the canonical 
and the microcanonical ensemble: 



(<^) CN - (^Ne) 



r 9/2 



MC 



/_^\i (k B TV< 

315C (|) \Hu) \ fruj J 



(70) 

It is easy to observe that the difference is proportional 
to N in the thermodynamic limit. Therefore the micro- 
canonical fluctuations, similarly to the canonical ones, 
are anomalous, and in the thermodynamic limit predic- 
tions of both ensembles become equivalent. 



V. CONCLUSIONS 

In this paper we have analyzed the fluctuations of the 
number of condensed particles in a weakly interacting ho- 



mogeneous and trapped Bose gases. For a homogeneous 
system we apply the particle-number-conserving formu- 
lation of the Bogoliubov theory, and obtain the analyti- 
cal results for the fluctuations in the canonical and mi- 
crocanonical ensemble, including the corrections due to 
the finite size of the system. Our derivation is based on 
contour integral representation in terms of spectral Zeta 
function. We determine the poles of Zeta functions re- 
sulting from the Bogoliubov spectrum, and develop their 
asymptotic expansions. This allows us to find analytical 
formulas for the mean number of noncondensed particles, 
its fluctuations and the mean energy, which are valid for 
arbitrary ratio of the temperature to the chemical po- 
tential. We determine the microcanonical fluctuations 
from the thermodynamic identity which relates them to 
the quantities calculated in the canonical ensemble. Our 
analysis shows that both the microcanonical and canon- 
ical fluctuations exhibit anomalous scaling, and in the 
thermodynamic limit fluctuations in these two ensem- 
bles become equal. A similar behavior is observed for a 
trapped gas. In this case we carry out our calculations 
within the standard Bogoliubov theory, and evaluate only 
the leading order behavior in the thermodynamic limit. 



Our results have been derived in the regime of temper- 
atures much lower than the critical temperature, where 
the effects of the quantum and thermal depletion on the 
population of the condensate can be neglected. We stress, 
however, that generalization of the obtained results for 
the case of higher temperatures, can be done in straight- 
forward way in the spirit of Popov theory [2g. To this 
end, for a homogeneous system parameter a should be 
replaced by ao = gno/A, with the density of condensed 
atoms no determined in a self-consistent way from the 
equation for the mean number of noncondensed parti- 
cles. In a similar way one can generalize the results for 
the trapped gas. In this case the total number of par- 
ticles N in equation for \x should be replaced by the 
number of condensed atoms Ao, with Nq calculated in 
a self-consistent a manner. 



Finally we would like to comment on the application of 
the particle-number-conserving method to calculate the 
fluctuations. For a homogeneous system we have verified 
that both conserving and nonconserving theories lead to 
the identical results. Such insensibility of the condensate 
statistics, can be attributed to the fact, that below the 
critical temperature, the condensate acts as a reservoir 
of particles for the noncondensed states. Thus, as long as 
the temperature is not too close to the critical tempera- 
ture, the probability distribution of the number of non- 
condensed atoms remains insensitive to the actual num- 
ber of atoms. The total number of atoms influences only 
the energy spectrum in the interacting gas. Therefore, 
predictions of conserving and nonconserving theories are 
the same with respect to the statistics of noncondensed 
part of the system. 
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APPENDIX A: RESIDUES OF SPECTRAL ZETA 

FUNCTION FOR THE BOGOLIUBOV 

SPECTRUM 

In this appendix we evaluate residues of the spectral 
Zeta function, which appears in the studies of fluctua- 
tions in a weakly interacting Bose gas: 



n#0 



-*/2 



(Al) 



The series in (|A1|) is convergent provided that Rez > |. 
In the remaining part of the complex plane, 9(z, a) is 
defined through analytic continuation. We rewrite Eq. 
(|A1|) in the following way 



the series 9>(z, a) containing rest of the terms 

9(z,a) = 9 < (z,a) + 9 > (z,a), (A6) 

9<(z,a) = J2 (n 4 + 2an 2 p /2 , 



0<n 2 <2a 

? > (z,a) = Y^ (n 4 + 2an : 

n 2 >2o 



2\-*/2 



(A7) 
(A8) 



The finite sum 6<(z,ot) is analytic on the whole com- 
plex z-plane. For the second term 9 > {z, a) we apply the 
binomial expansion, and obtain 

k=0 { 2 } 

x (Ce(z + k) - C, E {z + lb; 2a)) , (A9) 



where 



( E {z;a)= ^ (n 2 

0<n 2 <a 



(A10) 



Since the function Ce(z~, a) is analytic on the whole com- 
plex z-plane, the residues of 9{z, a) are given by Eq. 
1511. also in the case 2a > 1. 



9(z,a) = J2(n 2 



(1 + 2a/n 2 



3/2 



(A2) 



For the moment we assume that la < 1. Later, we show 
how the following derivation can be generalized for larger 
values of a. Since n 2 > 1, in the regime 2a < 1 we can 
apply the binomial expansion 



(1 + 2a/n 2 ) 



-2/2_^( -i) fc r(f + fc)( 2a ) fc 
*! r(f) {n 2f 



k=0 



(A3) 



Substituting this expansion into l|A2|l . and interchanging 
the order of summation, we arrive at the following series 
representation for 9(z, a) 



9(z, a) 



E *# E ^ (*')^-(--- >>■>■ 



fc=0 



k\ 



r(f) 



(A4) 



where C,e{z + k) denote Epstein Zeta function. Since the 
singularities of r(| + k) in the numerator occur exactly 
at the same points as the singularities of T(|) in the 
denominator, the only poles in the series (|A4|) . are the 
ones of (e(z + k). Epstein zeta function C,e{z) has only 
single pole at z = §, with residue equal to 2tt. Hence the 
poles of 9(z, a) are located at z = | — k for k = 0, 1, 2, . . ., 
with residues 



Res z 9(z,a) = 2^-^^ |{(2a) fe . 



fc! 



r(I 



(A5) 



In the case 2a > 1 , we can split the series I|A1|I into finite 
sum 9 < (z,a), containing terms for which n 2 < 2a, and 



APPENDIX B: ASYMPTOTIC EXPANSION OF 

SPECTRAL ZETA FUNCTION FOR THE 

BOGOLIUBOV SPECTRUM 

In this appendix we develop asymptotic expansion of 
spectral Zeta function 9(z,a). To this end we apply the 
method described in [33]. First, we rewrite IjAlfl in the 
following way 

9{z,a) = J2(n 2 y^ (n 2 + 2a)~ f , Rez > |. (Bl) 

Making use of the integral representation of Gamma 
function, Eq. (|B1|1 can be transformed into 

1 W n#0 Jo 

(B2) 
For the term e _tn we use integral representation of the 
Mellin-Barnes type 



ds_ 

2~ni 



T(s)(tnT s , 



(B3) 



with c € R and c > 0. Interchanging the order of in- 
tegration and performing the integral over t we arrive 
at 



9(z,a) 



r( 2 



1 ^ f c+l °° ds , N/ 2N 



s-z/2 



xr(|- J} )(2a)-" ' 2 . 

(B4) 
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By choosing c > 3/2— Rez/2, the sum over n is absolutely 
convergent and may be expressed in terms of Epstein zeta 
function 



1 rc+Joo ds z 

z ' a) = mU_^ r(s) H s + 2 



2/ J c— too 



xr(|- S )(2ar 2 / 2 . 



(B5) 

For large a the function under integral rapidly decay for 
Res < as long as s is not too large. When |Res| becomes 
comparable to a the integrand starts to grow, and finally 
becomes infinite for Res — ► — oo. Therefore, we close the 
contour of integration in 1B5|) on the left side, but include 
only a few leading order poles in the calculation of the 
integral. If we choose c < Rez/2, then the integrand has 
poles at s = 3/2 — z/2 (pole of the Epstein zeta function) 
and at s = — k for k = 0, 1,2, . . . (poles of T(s)). For 
z y^ 3 + 2n with n = 0,1,2,... all the poles are of the 
order of one and the asymptotic expansion reads 



8(z,a) -2vr 



£(§_- !)r(* 

M-l 



(2a) 



3/2- 



r(f) 
(-i) fc r(f + fc)CMj-fc) 

iq k\ r(f) (2a)*+*/3 

0(a- M ~ z/2 ), z^3 + 2n. 



E 



(B6) 



APPENDIX C: DETAILS OF CALCULATION 
FOR THE CANONICAL FLUCTUATIONS 

In this appendix we calculate the canonical fluctua- 
tions from the integral representation 

rc+ioo i 

<<5 2 iV e ) CN = / — Z(t, z, a)((z - l)T(z), 

+ S 2 N q , (CI) 

where spectral Zeta function Z(t, z, a) is given by 

Z(t, z, a) = t z (6{z, a) + 2a 2 6{z + 2, a)) , (C2) 

The poles of the integrand give the following contribu- 
tions: 

Pole of Q(z — 1) at z — 2. Its residue is equal to 1, 
and this pole gives the contribution Z(t, 2, a)T(2). Using 
asymptotic expansion of 6{z 1 a), we obtain 



Z(2 :.n in 2l-t a ( Jcb(2) + 



2V2a 



(C3) 



Pole of r(z) at z = 0. Its residue is equal to 1, and 
this pole gives the contribution Z(t, 0, a)C(— 1). Applying 
asymptotic expansion of 9{z, a) for large a we extract the 
leading order behavior 



Z(i,0,a)C(^l)^- I 2(2a) 3 / 2 . 



(C4) 



Poles ofT(z) atz= -l-2n, for n = 0, 1,2,.. .. Their 
residues are equal to (-l) 1+2n /(l + 2ra)!. Since C(-2 - 
2n) = 0, these poles do not give any contribution. 

Pole of r(z) at z = —2. Its residue is equal to 1/2 
and this pole gives the contribution Z(t, —2, a)((— 3)/2. 
With the help of the asymptotic expansion we find that 
this term is of the order of 0(N°), and therefore may be 
neglected in the final result. 

Poles of T(z) at z = -2 - 2n for n = 0,1,2,.... 
Their residues are equal to (— 1) 2+2 "/(2 + 2n)\, and the 
contribution from these poles is Z(t, —2 — 2ri, a)£(— 3 — 
2n)(-l)- 2 - 2 ™/(2 + 2n)\. However, Z(t, -2 - 2n, a) = 0, 
as can be deduced either from the binomial expansion 
(|A4(I or the asymptotic expansion l|B6|) . 

Poles of9(z, a) at z = 3/2 — n for n = 0, 1, 2, . . .. Their 
residues are given by l|A5(l . The contribution, which we 
denote by c\(t, a), from these poles is 



ci{t, a) 



71=0 



(-i)«r(| + f) 



m<*) n 



r(i-f) 

C(i-n)r(|-n). (C5) 

After some algebra we rewrite ()C5|) in the following way 



X 



Cl (t,a) = 2irt 3 / 2 J2- 



■})!(. 



n=0 



V27Tt/ 



V re C(n+I) 



xT(n+i)(i-n). (C6) 

Assuming a/2nt < 1, the summation can be done ana- 
lytically with the help of summation formula 33] 



Ch(o,, 1 -a;) 



^ k\ 

k=0 



x k T(a + k) 



((a + k), 



(C7) 



valid for |x| < 1. Hence, the final result is 



d(t,a) = {irtf /2 



2-KtJ 
a f3 
2^ H 



(H[- r l 



+ 



1 



- 

2vrf/ 



(C8) 



Since the function c±(t, a) remains analytic for a/2irt > 
1, validity of the result l|C8|l can be extended for a/2nt > 
1 by analytic continuation. 

Poles of6(z + 2,a) at z = -1/2-nforn = 0,1,2,.... 
Their residues are given by l|A5|l . The contribution 
Ci(t, a), from these poles is 



a 2 cx 



Wr;(-i) B r(| + §) (2a\ 



t 1 / 2 ^ 

n=0 



,1 



r( 



t 



2 

xCH-n)r(-I-n). 

The sum can be rewritten in the following way 

1 x ' ^ n C(- + |)r(n + |) 



(C9) 



c 2 (£, a) 



irt 1 / 2 ^ 

n=0 



a 

2-Ktj 



n\{n + \) 



(C10) 
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The summation can be simplified with the help of the 
following identity, which can be derived from l|C7f) 

^ x k r(| + k) C(| + k) 2 

gfc! r(|) fc + i = 3 5l(a;) ' N<1 ' (C11) 

where 



. . y^ 3fc + 2x 
«1W = 2^ fc27irZ^3/2 



fc=i 



fc 2 (fc + a;) 3 / 2 



(C12) 



Assuming a/2itt < 1, the final result for C2(i, a) reads 

.,2 



c 2 (i, a) 



2V^rt V27ri/ 



(C13) 



Since the function si(ir) remains analytic for x > 1, va- 
lidity of the result 1|C13|) can be extended for a/2iit > 1 
by analytic continuation. 

Finally we calculate the quantum fluctuations S 2 N q — 
2J2k^o u k v iL- Expressing S 2 N q in terms of 0(z,a) and 
applying asymptotic expansion l|B6(l . we obtain 



(C14) 



S 2 N q = ^9(2 7 a)o,^(2a)^ 



Summ ing up the contributions of Eqs. I|C3|) . (|C4|) . (|CJ8|> . 
(|C13|l . and l|C14(l we obtain the canonical fluctuations 
(GDI. 



APPENDIX D: OVERLAP INTEGRAL 
BETWEEN «(r) AMPLITUDES 

In this appendix we calculate the integral 
/ d 3 ru* nlm (r)u n w m '{r). For hw < n, we apply the 



approximation u n i m w f n i rn+ /2 where /„i m+ is defined 
by Eq. 1)63(1 . Performing the integration over the angular 
part of the function u n i m , we arrive at 



dr <im (*>«'/' 



^(4n + 2Z + 3)(4?i' + 2Z' + 3) 

X fimm'Sll'Inn'li (Dl) 



where 



/„„</ = / d yy 2 ' +2 (l- y 2 )Pr i/2 ' 0) (l-2y 2 ) 

xP^ +1/2 ' 0) (l-2y 2 ). (D2) 



To perform the latter integral, first we express 
P r [ l+1/2 '°\ X ) in terms of Pi l+l,2 ' l \x) and P^ M (x) 
utilizing recurrence relations for the Jacobi polynomials, 
and then perform the resulting integrals with the help of 
[34|. The final result reads 

T (n+l)(n + / + |) 

2(2n + Z + |) (2n + ; + f) 

I q " 77, 77 ' 

2(2n + / + |) (2n + / + i) 

(n + 1) (n + j + |) S n+1 , n , 
2 (2n + I + |) (2n + Z + |) (2n + I + f ) 

n (n+ I + \) S n -l,n' 
+ 2 (2n + / - |) (2n + / + |) (2n + j + |) 



(D3) 
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